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SUMMARY 


A numerical study was made to examine the effect of a porous surface on the 
aerodynamic performance of a transonic airfoil. The pressure jump across the normal 
shock wave on the upper surface of the airfoil was reduced by making the surface 
below the shock porous. The weakened shock is preceded by an oblique shock at the 
upstream end of the porous surface where air is blown out of the cavity. The lambda 
shock structure shown in the numerical result qualitatively agrees with that 
observed m the wind tunnel. According to the present analysis, the porous airfoil 
has a smaller drag and a higher lift than the solid airfoil. 


INTRODUCTION 


The transonic flow past an airfoil inflight at a high subsonic Mach number is 
characterized by a normal shock wave standing on the upper surface. If the airfoil 
surface in the neighborhood of the shock wave is made porous with a cavity under- 
neath it (as sketched in fig. 1), the pressure difference across the shock induces a 
flow through the cavity. The passive blowing and suction through the upper surface, 
before and after the shock respectively, reduces the strength of the shock. It can 
be expected that the wave drag of the porous airfoil is decreased. 

This expectation has been proven by the experimental results obtained by Bahi 
et al. and Nagamatsu et al. (refs. 1 and 2) which indicate that at high subsonic 
Mach numbers a properly arranged porosity can reduce both wave and viscous drag. In 
these experiments, however, the airfoil model was mounted so that the lower half of 
the airfoil was embedded in the bottom wall of the test section and only the upper 
surface was exposed to the transonic flow. This experimental condition is different 
from free flight. The negligible loss in lift on the porous airfoil reported in 
references 1 and 2 might not be the appropriate conclusion for a lifting airfoil. 

Numerical computations of the flow about a porous NACA 0012 airfoil were per- 
formed by Savu and Trifu (ref. 3) by solving the transonic small-perturbation poten- 
tial equation. Their results, for the airfoil at zero incidence and with the same 
porosity distribution on both the upper and lower surfaces, show that the standing 
shock may be eliminated completely by choosing a proper distribution of porosity for 
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a given Mach number. Additional references on thin airfoils with porous surfaces in 
pure subsonic or pure supersonic regimes can be found in reference 3. 

The aforementioned experimental and theoretical results indicate that porous 
surfaces can be used to reduce the wave drag of an airfoil at transonic speeds. 
Unlike a supercritical airfoil whose fixed shape is designed for certain optimum 
Mach number and angle-of-attack ranges, a porous airfoil can be adjusted without 
changing its contour shape. This adjustment is achieved with variable surface 
porosity and may result in improved performance over a wide range of conditions. 

The aerodynamic performance of a symmetrical airfoil in inviscid transonic flow 
with a partial porous surface is examined in the present work. Instead of the 
small-perturbation equation with linearized boundary conditions used in reference 3, 
the transonic full-potential equation is solved with boundary conditions applied 
exactly on the airfoil surface. The analysis also includes the influence of angle 
of attack which was ignored in previous computational studies. 

The work of the first two authors was supported by the U.S. Air Force Office of 
Scientific Research under Grant AFOSR 84-0037. Their use of the office and comput- 
ing facilities at NASA Ames Research Center is highly appreciated. 


ANALYSIS 


Because the entropy jump across the shock is third order, inviscid transonic 
flow is approximated as an isentropic and irrotational flow. This assumption allows 
the flow to be described by a potential $ which satisfies the continuity equation. 

( p 4 , x ) x + (P'i’y)y = 0 (la) 
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[- 


X_I_1 
Y + 1 




1/(y-1) 


(1b) 


where x and y are Cartesian coordinates nondimensionalized by the airfoil chord 
length c; p and <)> , ij> are the density and velocity components nondimensionalized 
by the stagnation $ens¥ty p and the critical sound speed a*, respectively; and 
Y is the ratio of specific fieats. Equation (la) (the full-potential equation) is 
written in strong conservation law form which is an important consideration for 
shock capturing algorithms. 

For computational convenience, these equations are transformed from the physi- 
cal domain (Cartesian coordinates) to a computational domain by the transformation 
(see fig. 2) 


C = C(x,y) n = n(x,y) 


(2) 
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The full-potential equation written in the computational domain (5-n coordi- 
nate system) is given by 


(pU/J) + (pV/J) n = o 


(3a) 
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The variables U and V are the contravanant velocity components along the 
j; and n directions, respectively; A^, A 2 , and A^ are metric quantities; and J is 
the Jacobian of the transformation. 

The body-fitted grids used m this work are generated by the finite difference 
solution of Poisson's equations using a computer code (GRAPE) developed by Sorenson 
(ref. 4) and based on the work of Steger and Sorenson (ref. 5). This grid genera- 
tion code allows control of both grid point spacing along and normal to the bound- 
aries and the angles at which grid lines intersect the boundaries. The C-mesh grid 
topology used in this study is shown in figure 2 . 

The governing equations (eqs. (3)) are solved using a fully implicit approxi- 
mate factorization scheme called AF2. A modified version of the TAIR computer 
program is used for this purpose (refs. 6 and 7). This code provides rapid conver- 
gence requiring only a few seconds of computer time per case on the CRAY XMP com- 
puter. The supersonic regions of flow are stabilized using an upwind bias of the 
density. This provides an efficient and reliable spatial differencing scheme for 
the capture of weak (transonic) shock waves. For more information about the flow 
solver algorithm or the TAIR computer program see references 6 and 7. 
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The boundaries associated with the physical domain are transformed to bound- 
aries of the computational domain. At the outer boundary, the velocity potential is 
set as the sum of the uniform free-stream component and the component due to a 
vortex with circulation r, where r is the jump in the velocity potential at the 
airfoil trailing edge. At the airfoil surface, there are two types of boundary 
conditions depending on the surface property. 

1. On solid regions of the airfoil, the flow tangency condition is satisfied 
by requiring that the contravariant velocity component in the n-direction vanish. 

V = 0 on solid surface (4) 

2. On porous regions of the airfoil, V does not vanish and its value is 
determined in the following manner. Through the porous surface the transpiration 
velocity v n is governed by Darcy's law such that 

v n = -(1/2)o(C p - C p ) (5a) 


oU /a* 


(5b) 


where the subscript n indicates the direction of the outward normal on the sur- 
face; Cp and C are the pressure coefficients above and below the porous surface, 
respectively; o is the porosity distribution function which is determined by vis- 
cosity as well as by the size and density of the holes in the porous surface. For 
passive flow through the cavity, the net mass flow through the porous surface of 
length s must be zero, or 


- = / » v n ds = 0 

¥ s 


which gives, by assuming a constant cavity pressure, 


C = I opC ds / I op ds 

P J s P / J s 


( 6 ) 


The constant cavity pressure assumption was made for convenience but is expected to 
be qualitatively accurate. In addition, several variations of porosity have been 
investigated with the expectation that this would simulate a variable cavity pres- 
sure. It is planned that additional work to improve the cavity flow numerical model 
will be completed in the near future. 

The transpiration velocity v n computed from equations (5) is then transformed 
into the computational domain (ref. 8): 
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With the C-mesh topology, the n-coordinate lines intersect the body at close to 
right angles, and to a good approximation = 0 at the body. Using this approxi- 
mation v n does not contribute to U and can be expressed in the computational 
domain as: 


V = A 3 v n on the porous surface (7b) 

The numerical implementation of the boundary conditions on the airfoil surface 
is outlined as follows: 

1. On the solid portion of the airfoil, the boundary condition of flow tan- 
gency (i.e., V = 0) is used to obtain 

(r) = -if) (8) 

J l ,NJ+( 1/2) J l ,NJ-{ 1/2) 

where j = WJ is the airfoil surface (see fig. 2). 

2. On the porous portion of the airfoil, the boundary condition (i.e., 

V = /^ 3 v n ) is applied with the aid of Taylor series expansions 


(oV\ _ /pV\ + An / pV\' (An) 2 (Pi 

J i,NJ±( 1/2) ' J l ,NJ 2 l ,NJ 8 J 


II 

i,NJ 


± 0[ (An) 3 ] . 


so that the first-order boundary condition becomes 


(9) 



i,NJ+( 1/2) 



i,NJ-( 1/2) 



( 10 ) 


In fact, equation (8) is a special case of equation (10) for V = 0 at j = NJ. 

For more information on the transpiration velocity boundary condition see refer- 
ence 8. 

Equations (3), together with the surface boundary conditions (eqs. (8) 
and (10)) and appropriate outer boundary conditions, are solved numerically using 
the modified TAIR computer code. At the beginning of each iteration, the porous 
boundary condition is updated using equations (5). 

Once the flow field is obtained, forces on the airfoil in the x and y direc- 
tions can be computed by performing integrations around the airfoil contour. 
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( 11 ) 
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( 12 ) 


in which S is the circumference of the airfoil; u and v are velocity components 
in the x and y directions, respectively; and n x and ny are direction cosines 
between the airfoil surface and the x and y directions, respectively. Equa- 
tions (11) and (12) are derived from the momentum conservation theorem. Note that 
the second term in both equations (P v n u and pv n v) are zero for solid wall airfoils 
but must be retained in the present calculations because of the porous wall assump- 
tion. Once F x and F are obtained from equations (11) and (12), the lift and drag 

coefficients, and Cp, are easily computed. 


NUMERICAL EXPERIMENTS AND RESULTS 


A NACA 0012 airfoil was used for a series of numerical studies to investigate 
the effect of a porous surface on its transonic flight performance. In all of the 
cases presented here, computations were carried out on a 223 * 31 C-type mesh, 
having 162 of the 6913 grid points placed on the airfoil surface. 

Three types of porosity distribution have been examined which were obtained by 
varying the porosity distribution function o in equation (5b). They are described 
as follows: 

Type 1 

a = 

Type 2 


a = 


the Xi and ^ parameters are the limits of the porous region shown m figure 1, 
and ° max 1S the maximum porosity at the midpoint of that region. The Type 2 
porosity distribution is the same distribution function used by Savu and Trifu 
(ref. 3). 


constant 


max 


^ sin 


x - x. 


(13) 


(14) 


6 



Type 3 



The x„ parameter is the horizontal position of the shock that would result if the 
porous surface were solid, and x k represents either x^ or Xg, depending on 
whether x is less or greater than x g . This function automatically adjusts the 
porosity distribution so that it decreases from the maximum value o max under the 
shock wave to zero at either end of the porous region. Numerical experiments indi- 
cate that, in most cases, the Type 2 and 3 distribution functions are superior to 
the uniform a distribution for drag reduction. For this reason, the results 
obtained based on the Type 1 distribution function are not reported here. 

The effect of a porous surface on a Mach 0.8 flow past the NACA 0012 airfoil at 
zero angle of attack is presented in figures 3 and 4. A Type 2 porosity is distrib- 
uted on the upper surface between = 0.4 and x 5 = 0.8 with a =0.6. The 
constant-Mach number contours show a normal shock standing below the solid lower 
surface of the porous airfoil whose location and strength are not much different 
from those of the shock on the original solid airfoil. However, the shock above the 
upper surface is weakened in the presence of the porous surface, as shown by a group 
of less concentrated constant-Mach lines around the sonic line. The figure also 
reveals that the porous surface causes a weak oblique compression wave at its 
upstream end ( x -j = 0.4), due to the blowing from the cavity, and a readjusted com- 
pression downstream of the shock due to the suction of air into the cavity. The 
contour lines in the shock region are no longer normal to the airfoil surface, 
resulting in a lambda-shaped shock wave structure similar to that photographed m 
the laboratory (refs. 1 and 2). 

The pressure distribution on the upper airfoil surface for the case just pre- 
sented is plotted in figure 4. The results for both a porous airfoil (dashed line) 
and solid airfoil (solid line) are compared. The comparison clearly shows that the 
original steep compression through the normal shock on the solid airfoil has been 
reduced in the presence of the porous region. The original shock is replaced by 
several weaker compressions over the region covered by the porous surface. The 
resulting weaker adverse pressure gradient would lessen the possibility of flow 
separation. The porous upper surface has a negligible influence on the pressure 
distribution along the lower surface; that distribution is, therefore, omitted 
here. The asymmetrical pressure distribution on the upper and lower surfaces causes 
a small lift on the airfoil at zero angle of attack, with C L = 0.0183. On the 
other hand, a 27 . 5% decrease in wave drag from C D = 0.0069 to 0.0050 is found for 
the porous airfoil. 

Much larger decreases in drag are obtained when both upper and lower surfaces 
are made porous. The effect of varying the porosity strength on the drag of a 
double-porous NACA 0012 airfoil at a = 0° is shown in figure 5. A Type 2 porosity 
distribution between x^ =0.1 and X£ = 1.0 is used for all curves presented. 

The porous surface has a drag-reduction effect only when a shock appears above it at 
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Mach numbers higher than 0.77, and that effect is enhanced by increasing o 
(i.e., by using porous surfaces having smaller resistance to the penetrating a £low) . 

With the same a m , changing from Type 2 to Type 3 porosity does not cause a 
significant change in m lrag if the normal shock appears near the center of the porous 
surface; otherwise Type 3 porosity is considerably better than Type 2 in smearing 
the shock and reducing the drag. 

It is not possible to make a direct comparison of the computed data shown in 
figure 5 with the experimental results presented in reference 2 obtained in the wind 
tunnel under a different flow arrangement. Furthermore, the measured drag data 
contain the effects of viscosity, airfoil surface roughness, and other factors that 
have not been considered in our analysis. Nevertheless, the drag reduction effect 
caused by the porous surface observed in the laboratory is qualitatively the same as 
that obtained in the present study. 

The behavior of porous airfoils at angle of attack is now described. The Mach 
number contours for a solid NACA 0012 airfoil at Mach 0.75 and an angle of attack of 
1° are plotted in figure 6. A shock wave appears only on the upper surface. A 
Type 3 porosity of strength ° max =0.3 is then distributed on 90$ of the upper 
surface between Xi =0.1 and x 2 = 1.0. The resulting flow pattern and pressure 
distribution, plotted in figures 7 and 8, respectively, reveal that this widely 
distributed porosity is very effective in reducing the shock strength. By making 
the upper surface porous, increases from 0.2397 to 0.3566 while Cp decreases 
from 0.0024 to 0.0008, which corresponds to a nearly shock-free condition. 

To study the effect of varying porosity strength on lift and drag of a tran- 
sonic airfoil, a Type 3 porosity is distributed between x-| = 0.3 and x 2 = 0.9 on 
the upper surface of a NACA 0012 airfoil. For this series of cases the angle of 
attack is fixed at 1°. The result for drag plotted m figure 9 is analogous to that 
shown in figure 5 for the double-porous airfoil at zero angle of attack, showing 
that large drag reductions can be achieved by increasing the porosity. Plotted m 
figure 10 are the lift coefficient data which indicate that lift is increased by 
making the upper surface porous. Unlike the drag coefficient, is affected by 

the porous surface at Mach numbers less than 0.72 when the shock is still upstream 
of the porous region. The higher lift is caused by asymmetric changes in the pres- 
sure on the porous upper and solid lower surfaces of the airfoil. 


CONCLUDING REMARKS 


The numerical solution to the full-potential equation shows that the drag of a 
transonic airfoil can be reduced when a porous surface is placed in the shock 
region, and that the drag reduction ability increases with increasing porosity 
°max’ ^ or an air ^ 01 ^ having porous upper and lower surfaces at zero incidence, the 
results reported m reference 3 agree qualitatively with those obtained here without 
using the small-disturbance approximation. Computed flow fields display the lambda 
shock structure that has been observed in the laboratory. It is also found that the 


8 



porous upper surface alone causes not only a reduced drag, but also an increased 
lift on a transonic airfoil at angle of attack. 

The outcome of our preliminary study of using a porous surface for drag reduc- 
tion is encouraging. In our continued study we will refine the present potential 
flow model by computing the internal cavity flow instead of assuming that the cavity 
pressure is constant. Thus the shape of the cavity will have an influence on the 
drag, which is a fact already being found experimentally (ref. 1). Viscous and 
three-dimensional effects will also be implemented to simulate the flow more realis- 
tically. Finally, this work will be extended to handle both passive and forced 
suction and blowing on the airfoil surface. 


9 



REFERENCES 


1. Bahi, L.; Ross, J. M. ; and Nagamatsu, H. T.: Passive Shock Wave/Boundary Layer 

Control for Transonic Airfoil Drag Reduction. AIAA Paper 83-0137, Jan. 1983 

2. Nagamatsu, H. T.; Orozco, R. D.; and Ling, D. C.: Porosity Effect on Supercnt 

ical Airfoil Drag Reduction by Shock Wave/Boundary Layer Control. AIAA 
Paper 84-1682, June 1984. 

3. Savu, G.; and Trifu, 0.: Porous Airfoils in Transonic Flow. AIAA Journal, 

vol . 22, July 1984, pp. 989-991. 

4. Sorenson, R. L. : A Computer Program to Generate Two-Dimensional Grids About 

Airfoils and Other Shapes by the Use of Poisson's Equation. NASA TM-81198, 
May 1980. 

5. Steger, J. L.; and Sorenson, R. L.: Automatic Mesh-Point Clustering Near a 

Boundary in Grid Generation with Elliptic Partial Differential Equations. 
Journal of Computational Physics, vol. 33, Dec. 1979, pp. 405-410. 

6. Holst, T. L.: Implicit Algorithm for the Conservative Transonic Full-Potential 

Equation Using an Arbitrary Mesh. AIAA Journal, vol. 17, Oct. 1979, 
pp. 1038-1045. 

7. Dougherty, F. C. ; Holst, T. L.; Gundy, K. L. ; and Thomas, S. D.: TAIR-A Tran- 

sonic Airfoil Analysis Computer Code. NASA TM-81296, May 1 98 1 . 

8. Van Dalsem, W. R. : Simulation of Separated Transonic Airfoil Flow by Finite- 

Difference Viscous-Inviscid Interaction. Ph.D. dissertation, Stanford 
University, June 1984. 


10 




Figure 1.- Porous airfoil in transonic flow 
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Figure 3.- Constant Mach number contours around a porous NACA 0012 airfoil, 
= 0.8, a = 0°, Type 2 porosity, a max = 0.6, = 0.4, ^ = 0.8 (upper 

surface is porous). 
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Figure 5.- Effect of varying porosity strength on the drag of a NACA 0012 airfoil 
at a = 0° , Type 2 porosity, x^ =0.1, *2 - 1-0 (both the upper and lower 
surfaces are porous). 
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Figure 6.- Constant Mach number contours around a solid NACA 0012 airfoil for 

M = 0.75 and a = 1°. 
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Figure 7.- Constant Mach number contours around a porous NACA 0012 airfoil, 

M = 0.75, a = 1°, Type 3 porosity, o = 0.3, Xi =0.1, x ? = 1.0 (upper surface 
°° max ' ^ 

is porous). 




x/c 

Figure 8.- Pressure distribution around the airfoil described in figure 7, porous 
airfoil — dashed line, solid airfoil — solid line. 
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Figure 9 
at i 



- Effect of varying porosity strength on the drag of a NACA 0012 airfoil 
= 1°, Type 3 porosity, = 0.3, x 2 = 0.9 (upper surface is porous). 
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